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Abstract. Nanometer- and micrometer-sized solid particles play an important role in the evolu- 
tionary cycle of stars and interstellar matter. The optical properties of cosmic grains determine 
the interaction of the radiation field with the solids, thereby regulating the temperature structure 
and spectral appearance of dusty regions. Radiation pressure on dust grains and their collisions 
with the gas atoms and molecules can drive powerful winds. The analysis of observed spectral 
features, especially in the infrared wavelength range, provides important information on grain 
size, composition and structure as well as temperature and spatial distribution of the material. 

The relevant optical data for interstellar, circumstellar, and protoplanetary grains can be ob- 
tained by measurements on cosmic dust analogs in the laboratory or can be calculated from grain 
models based on optical constants. Both approaches have made progress in the last years, trig- 
gered by the need to interpret increasingly detailed high-quality astronomical observations. The 
statistical theoretical approach, spectroscopic experiments at variable temperature and absorp- 
tion spectroscopy of aerosol particulates play an important role for the successful application of 
the data in dust astrophysics. 
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1 INTRODUCTION 

Dust grains exist in a wide variety of cosmic environments ranging from the diffuse interstellar 
medium and molecular clouds in our galaxy to circumstellar shells and disks around young 
stellar objects, dusty outflows from evolved stars, atmospheres of brown dwarfs, the dust tori 
of active galactic nuclei and starburst galaxies, and high-redshift quasars. The particles are 
especially important during the birth process of stars and as precursor material for the formation 
of planets. In the late stages of stellar evolution they condense in molecular outflows and are 
important for the mass loss process and the spectral appearance of AGB and post-AGB stars. 

The field of cosmic-dust research has been the subject of intense research over the last 
decades. The international conferences "Astrophysics of Dust" in Estes Park, Colorado (2003) 
and "Cosmic Dust - Near and Far" in Heidelberg, Germany (2008) both reflected the rapid 
development of the field. Astronomical dust research is of course triggered by the progress in 
astronomical observations, mainly by those in the infrared range, because of the opportunities to 
(1) observe dusty objects optically thick in the visible and (2) to analyse the thermal emission 
of dust by infrared spectroscopy. The ISO and Spitzer space telescopes have played a major 
role in this respect, and current and future projects such as the Herschel mission and JWST will 
continue to ensure progress in this field. 

The description of the interaction of radiation fields with a system of solid particles always 
requires knowledge of their absorption and scattering cross sections which depend on their 
chemical composition, solid-state structure and morphology (see Table [T]). Any continuum 
radiative transfer code for the calculation of intensity/polarization maps and spectral energy 
distributions needs the specification of these properties. In addition, the transfer of momentum 



Table 1. Chemical and structural properties of solid particles and the description of optical 
properties on different length scales 



Length scale Structural properties 



Description 



"Atomic" 



"Mesoscopic" 



Chemical composition and bonding 
Crystal structure 
Defects and impurities 

Inhomogeneitie s 
Porosity 
Mantles 
Surface states 



Dielectric function e(A) 
Complex refractive index m(A) 



Theory of effective media 

ecff(A), 

Core-mantle description 



"Macroscopic" 
Morphology 



Size distribution 
Grain shapes 
Agglomeration 
Coalescence 



Mie theory, T-matrix method. 
Separation of variables method, 
Discrete dipole approximation 

C'ext(A), Cabs(A), Csca('^) 



from a radiation field to a grain is determined by the extinction and scattering efficiency and the 
average cosine of the scattering angle. 

In this review we will deal with the progress in theoretical and experimental work aimed to 
providing data and tools for the interpretation of spectroscopic dust-related observations, and 
for modeling of dusty astrophysical objects. Current models go significantly beyond the as- 
sumption of spherical grains characterized by artificially constructed optical constants, as has 
been the prefered approach for many years. Effects of grain morphology, temperature, and 
composition have been studied experimentally, and complex light-scattering models have been 
used to explore the dependence of scattering and absorption cross sections on the structure of 
grains on all length scales (see Table [U, including porosity, inhomogeneity, and agglomera- 
tion. Such grain models have been used in sophisticated multi-dimensional radiative transfer 
codes. Among the experimental data, which can be included as inputs into the computations, 
are also measurements of the scattering matrices of cosmic dust analogs which have been per- 
formed in Amsterdam for many years and have been made available through the Amsterdam 
Light Scattering Database ( http://www.astro.uva.nl/scatter) l. Although this review unfortunately 
cannot cover this and some other experimental fields devoted to interpret observational data of 
cosmic dust, all general considerations presented in the next sections are valid for them as well, 
especially those on the relation between calculations and measurements. 

It is beyond the scope of this paper to discuss extensively the astronomical dust models and 
optical data derived from observations. For a description of the various "classical" dust models, 
we refer to the contributions in Ref. 1. More recent models have been developed by Mathis fj], 
Li & Greenberg Ul, Zubko et al. 21, Draine & Li fS) and can also be found in the proceedings 
of the two dust conferences ||6]|7l . 

2 TWO WAYS FOR PREDICTING OPTICAL DATA OF COSMIC DUST 
PARTICLES 

Our knowledge about the nature of cosmic dust grains mainly comes from spectroscopic obser- 
vations over the infrared, optical, and UV wavelength ranges. Infrared spectroscopy provides 
information about vibrational modes of functional groups and crystal lattices, whereas the opti- 
cal and UV range of the galactic extinction curve is mainly dominated by a strong UV resonance 
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Fig. 1. Two ways for determining optical properties of cosmic dust analogs. 

at 217 nm, usually attributed to an electronic transition in carbonaceous grains or poly aromatic 
hydrocarbons (PAHs). In the infrared, features related to the abundant silicates are frequently 
observed, especially the bands around 10 and 18 fim attributed to Si-O stretching and bending 
vibrations. 

Other constraints for the chemical and physical structure of grains are provided by abun- 
dance considerations, by our understanding of the temperature and density range, as well as 
by the presence of shocks and winds in interstellar and circumstellar media (see, e.g., Refs. 8 
and 9). In oxygen-rich environments we can expect the presence of oxide materials (mainly 
silicates) and oxygen-containing ices (H2O, CO2). Carbon-rich environments contain mainly 
carbonceous solids, carbides, and sulfides. Here, we will only deal with refractory materials 
and not with molecular ices lllOlll II and larger molecules such as PAHs 112111311 . 

After the discovery of abundant crystalline minerals in circumstellar disks (e.g., Refs. 14 
and 15), and the increasing efforts to study the detailed mineralogy and morphology of the 
dust in these environments from its infrared emission spectra, both theorists working on light- 
scattering theory and experimentalists in astrophysical laboratories have focused their work on 
predicting the spectroscopic properties of cosmic dust particles more accurately. In contrast to 
the amorphous silicates in the diffuse interstellar medium, crystalline silicates are frequently 
observed in protoplanetary disks, which indicates their in-situ formation by thermal annealing 
or shocks (for a review see Ref. 16). An interesting result of astronomical observations of pro- 
toplanetary and circumstellar dust is the fact that the crystalline silicates are mostly magnesium- 
rich IlTHIll. 

Extracting reliable dust properties from the thermal emission spectra of a mixture of dust 
particles requires precise knowledge of the optical properties for all possible constituents under 
various conditions. Crystalline silicates as the most important dust minerals, e.g., have very 
characteristic infrared spectra with many bands at wavelengths between 8 fim and the submil- 
limeter range. The details of band positions and band profiles, however, depend on the chemical 
composition, the grain morphology, and the grain temperature. This means that morphological 



information, e.g., can only been extracted from the spectra if the influence of the other parame- 
ters is sufficiently well understood and can be disentangled. 

The grains condensed in circumstellar environments and modified by interstellar processes 
and even re-formed in the interstellar medium may be quite chaotic in composition and struc- 
ture. Table[T]summarizes some of the chemical and structural properties which may characterize 
the particles on different length scales. The theoretical description of the optical properties of 
such particles can be very difficult. Even the "optical constants" (the complex refractive index 
m = n + ik OT the dielectric function e = ei + 162) of the materials constituting the particle 
may be poorly known. 

These facts require studies which have to be done in the laboratory, via synthesis and optical 
spectroscopy of materials representing important aspects of the cosmic solid matter 1201 12TI . 
The investigations of these "cosmic dust analogs" do not necessarily aim at exactly reproducing 
cosmic dust. In our terminology dust analogs comprise also bulk matter or thin films prepared 
for the derivation of optical material constants. In addition, we also consider simple "model" 
particles (e.g., silica monospheres) used for the investigation of special aspects of particle optics 
(e.g., agglomeration effects). 

Figure[T]illustrates the various aspects of this approach. Two principal possibilities exist for 
the derivation of the optical properties of particles: 

1 . Calculation of cross sections from optical material data which may have been measured 
or taken from a data base. In most cases, the data will originate from a bulk or thin-film 
material since measurements on these geometrically simple systems are much easier and 
much more work has been done on them. Therefore one has the problem that the material 
may not be exactly the one one would like to have. Special material structures may only 
exist in the form of small grains (e.g., silicon nanoparticles, carbon onions). Additionally, 
theory has some problems with the calculation of spectra for complicated grain shapes in 
strong absorption bands. This will be discussed further in Sec. 4. 

2. Direct measurement of the optical properties on particle samples. This certainly is the 
preferred but also experimentally much more demanding possibility. One has to pro- 
duce particles of the desired composition and structure and one has to store them in an 
environment as similar to the cosmic one as possible. This is sometimes problematic be- 
cause conventional preparation of particulate samples for spectroscopy (i) uses matrices 
in which the particles are embedded and (ii) is not able to prevent agglomeration of the 
primary particles to larger clusters. 

The next section will treat the basics of the first (computational) way, the light scattering 
theory, which also underlies the experimental way. Section 4 will deal with measurements of the 
material optical data from bulk and thin film samples as well as of extinction data of particulate 
samples. 

3 ASTRONOMICAL FEATURES AND DEVELOPMENTS IN LIGHT 
SCATTERING THEORY 

The optical properties of small particles can considerably deviate from those of bulk materials 
because of the occurrence of surface modes. The structure of the interface of small particles 
with the surrounding medium usually has strong effects on their optical properties. In addition, 
the shape and internal structure of the particles are important parameters in determining the 
position and shape of the observed astronomical features. Furthermore, the optical depth and 
the temperature structure of the dusty astronomical media are decisive for the occurrence of 
the solid-state bands and will determine if they are observed in emission, absorption, or self- 
absorption. Protostars, deeply embedded in their parental molecular clouds, show absorption 



features due to molecular ices (e.g., H2O, CO, CO2 and NH3) and refractory silicates. In 
contrast, the spectra of protoplanetary disks around young stars are mostly characterized by a 
variety of silicate emission bands, produced in the optically thin disk atmosphere on top of the 
optically thick disk region. 

For micron-sized grains classical electrodynamic theory can be applied using bulk optical 
constants. Textbooks which treat the classical electrodynamics of absorption and scattering of 
light by small particles are the publications by van de Hulst 1221 . Kerker 1231 , Bohren & Huff- 
man Il24l . Kokhanovsky ||25l, Mishchenko, Hovenier & Travis I^S], and Mishchenko, Travis & 
Lacis |l27l . 

The qualitative features of the absorption and scattering of light strongly depend on the ratio 
between the wavelength of the incident light A and the typical radius of the particle a. We can 
distinguish three regions: 

1 . A ^ a: This is the limit of geometrical optics, where the propagation of light is described 
by rays which are reflected and refracted at the surface of the scatterer according to Snell's 
law and the Fresnel formulae. For absorbing materials, light can penetrate only within 
the skin depth. Scattering, therefore, is mainly a surface effect and the extinction cross- 
section per volume goes roughly as 1/a. 

2. X ^ a: This is the regime of wave optics, where both the angular and the wavelength 
dependence of the cross sections are dominated by interferences and resonances. For 
spherical particles, these can be calculated by the Lorenz-Mie theory ll28l . 

3. A » a: This is the Rayleigh limit. If furthermore we have A » \m\a, where m is the 
(complex) refractive index of the particle, we are in the quasi-static limit. Then, both the 
incident and the internal field can be regarded as static fields. In this regime, phase shifts 
over the particle size are negligible. This implies that it is generally sufficient to consider 
only the dipolar electric mode. 

The interaction of infrared radiation with sub-micron sized cosmic grains can generally be 
considered as a good example of the quasi-static case. However, particles with high imaginary 
parts of the refractive index (metals, semiconductors) may easily violate this condition even 
at infrared wavelengths. The same is true for laboratory measurements where it is not always 
guaranteed that the particles fulfill the condition for the quasi-static case. 

For the sake of simplicity and physical insight, we will deal in the following only with the 
quasi-static case. The extinction cross section Ccxt is given by: 

Ccxt ~ Cabs + Cscaj (1) 

where Cabs and Csca are the absorption and scattering cross sections. According to Rayleigh's 
law Csca fc^ with the wavenumber k = 2tt/X. For A S> a this means in practice that the 
extinction cross section equals the absorption cross section. The cross section for the radiation 
pressure can be computed via the formula 

Cpic = Ccxt - Csca(cOs6'), (2) 

with {cos6) the mean of the cosine of the scattering angle. Following the conventions in ll24l . 
the expression for the extinction cross section is given by: 

Ccxt = fclm(a), (3) 

with a being the polarizability. This quantity depends on the complex dielectric function, e = 
ei + i 62, of the particle material and on the dielectric function of the embedding medium, em, 
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Fig. 2. Shape effects for Si02 particles in the quasi-static limit. CDE refers to a continuous 
distribution of ellipsoidal shapes (Eq. 9). 

which is in most cases a real wavelength-dependent number In addition, the polarizability is 
determined by the shape of the particle. 

For ellipsoids the polarizability, ai, in an electric field parallel to one of its principal axes is 
given by 

em + Li[e - eni) 

where Lj are geometrical factors and V is the volume of the ellipsoid. The relation Li+L2+Ls=l 
implies that only two of these factors are independent. In the case of prolate spheroids with 
equal minor axes of length b=c and with eccentricity e = -^/l — c^/a^ (a the length of the 
major axis) the geometrical factors are given by 

1 - / 1 1 + e\ 
Li = -1+ In- (5) 

e \ 2e 1 — e / 

and L2=i3=(l — Li)/2. In the other case of oblate spheroids we have a=b and Li=L2=(l — 
L3)/2 with 

The simplest case of light scattering in the quasi-static limit is that for a sphere where 
a = b = c and Li = 1/3. This gives ai=a 

a = 3V^^^. (7) 

These equations demonstrate that resonances for particles surrounded by a non-aborbing medium 
occur close to wavelengths where the imaginary part of the dielectric function is close to zero 
and the real part fulfills the condition 



(8) 



This immediately implies that the wavelengths of resonances depend on the shape of the parti- 
cles. The resonances can only occur in regions where ei is negative (for a sphere ei — ~2eni). 
Examples of astronomically relevant materials which have negative values of ei are SiC 1291 
and Si02 lf30l in the infrared and graphite in the UV ifSTl . For metal particles resonances even 
occur at wavelengths where the bulk material does not show any absorption band. In contrast, 
the lattice features of amorphous silicates do not show this behaviour. 

In the case of lattice modes, the resonances are located between the transverse and logi- 
tudinal optical phonon frequencies. As an example, we show shape effects for Si02 particles 
in Fig. 121 The equations also show that the positions of resonances depend on the surround- 
ing medium - a fact which is often considered in astronomy in the case of dust models with 
core-mantle grains. For particles of arbitrary shape there are no analytical solutions of the light 
scattering problem (even in the quasi-static limit) available. Numerical models frequently used 
for non-spherical particles are the separation of variables method, the T-matrix method and the 
discrete dipole (or multipole) approximation (see, e.g., Ref. 32). The treatment of anisotropic 
materials gets more complicated, but is still possible 1331 . 

In reahty, we always expect a distribution of particle shapes. The difficulty both in astro- 
nomical applications and in laboratory experiments is to evaluate this shape distribution. A 
simple model for a distribution of particle shapes is the continuous distribution of ellipsoids 
(CDE), which by averaging Eq. (4) over all L; leads to 

2e 

aav = V -Lege, (9) 

e — 1 

where Log e denotes the principal value of the logarithm of the complex number e (see Bohren 
& Huffman 1983). The CDE assumes equal probability for the presence of every shape and 
averages over all orientations. This model can be used to estimate in an approximate way how 
important shape effects for a special material could be. 

Recently, Min et al. 13411 have shown that an arbitrary grain shape can be described by a 
distribution of geometrical factors P(L), which can be derived from a spatial discretization of 
the particle shape. Such distributions containing the precise shape information can be averaged 
for representative shapes of an ensemble of particles and can be used to calculate the ensemble 
polarizability 

aav=V ( P{L)a{L)dL (10) 
Ja 

with q;(L) according to Eq. (4). This "distribution of form factors (DFF)" model has been 
proven to be able to successfully describe experimental IR spectra for different grain shapes of 
submicron-sized particles and to predict the influence of an embedding medium correctly ll35l . 
P(L) distributions - DFFs - can be derived by fitting experimental spectra and typical DFFs 
for irregular and roundish grains have been extracted (see also Sec. 4.3.). Compared to other 
models, the use of these DFFs allows a fast computation of quite realistic band profiles for real 
particulates. 

Here we should stress that all of the previous considerations are only valid in the quasi-static 
case and for isolated particles. Another shortcoming of this and all similar models is that for 
crystallographic anisotropic materials, the polarizability has to be calculated with the different 
sets of optical constants separately. This neglects possible interaction of dipoles with different 
orientations and can lead to errors in some cases ll35l . 

In dense regions of interstellar matter (dense molecular cloud cores, circumstellar disks and 
envelopes), the submicron-sized grains can coagulate and finally form larger particles II36II37I . 
For a description of the interaction of electromagnetic radiation with fluffy aggregates composed 
of individual particles two distinct approaches are possible: 



1 . "Deterministic" approach: The frequency-domain Maxwell equation is solved for an indi- 
vidual cluster The resulting cross sections are calculated for many clusters and averaged 
over the ensemble and the orientation of the clusters. An advantage of this approach is 
that for special systems (e.g., clusters of spheres), exact solutions of the problem do exist. 
For computational reasons, however, these methods are limited to either comparatively 
small clusters or moderately absorbing systems. For larger clusters of highly absorbing 
matter (as in the case of graphitic carbonaceous grains and metal clusters), one still has 
to rely on approximations. Examples of this kind of approach are the discrete dipole 
and multipole approximations (DDA/DMA) and the extended Mie theory for multisphere 
aggregates. 

2. "Statistical" approach: The equations are formulated in terms of statistically relevant 
quantities only (e.g., average radial density function of the clusters, density correlation 
function) without any expUcit treatment of the properties of individual particles. Whereas 
a given cluster generally does not have any symmetry, statistical averages show rotational 
invariance, unless alignment mechanisms (magnetic fields, winds) break this symmetry. 
The advantage of this approach is that only the necessary information (ensemble and 
orientation averaged quantities) enter the calculations. Examples of this approach are 
the different effective medium theories, the strong permittivity fluctuation theory, and the 
"Distribution of Hollow Spheres" model 1381 . The above-mentioned DFF model belongs 
also to this class. 

Benchmark results for these methods and astronomically relevant materials have been pub- 
Ushed by Stognienko et al. |[39|, Michel et al. EOl, Xu & Gustafson llH, and Voshchinnikov et 
al. 1321 . The interested reader should consult these papers for more detailed information. The 
properties of porous and fractal cosmic grains have recently found renewed interest and were 
calculated by various authors l42H45]| . 

The 10 /im silicate feature has been frequently used to trace the particle size in protoplane- 
tary disks. Infrared spectroscopy of disks have shown the presence of flat-top profiles, indicative 
of the presence of micron-sized particles which are larger than the typical particle sizes in the 
diffuse interstellar medium (e.g., Refs. 46 and 47). The influence of porosity and aggregate 
structure on this important astronomical dust feature has been investigated by Min et al. l48l 
and Voshchinnikov & Henning l49l . 

4 DEVELOPMENTS IN SPECTROSCOPIC EXPERIMENTS 
4.1 Spectroscopic data of dust minerals 

One type of information which may be possible to extract from infrared emission spectra of cos- 
mic dusty media is the chemical composition and mineralogical structure of the dust minerals at 
the present stage of evolution of the respective object. In a great effort to provide experimental 
dust spectra for comparison with the astronomical infrared data, the spectroscopic properties of 
the (likely) major dust minerals have been investigated in different laboratories and have been 
made available to the astronomical community. For example, the infrared spectra of crystalline 
silicates with different iron to magnesium ratios have been studied in great detail in the last 
decade, especially by the group of C. Koike at Kyoto Pharmaceutical University (e.g., Refs. 50 
and 51) by manufacturing and measuring synthetic crystals of variable composition. 

The technique applied for the infrared spectroscopy is usually transmission measurements 
of powders embedded into solid matrices (pellets made of KBr, Csl or polyethylene depending 
on the wavelength range of the spectroscopic measurement). In order to reach the Rayleigh 
case (size smaller than wavelength), the particles are mechanically ground and size-separated 
by sedimentation. The fine-grained particles are dispersed among the matrix material as homo- 
geneously as possible at a certain mass ratio. The pellets then contain the particles at a known 
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Fig. 3. Mass absorption coefficient of talc (Mg3[Si40io|(OH)2]) at far-infrared wavelengths for 
two different temperatures. The ratio 1:10 (by mass) refers to the dilution of the talc powder by 
polyethylene in the pellet — after Ref. 55. 

column density a = m/A, with A being the pellet cross sectional area. The advantages of this 
technique are easy preparation, small sample consumption and quantitative determination of 
mass absorption coefficients k = —ln{T)/a from the spectral transmission of the pellet. Due 
to Kirchhoffs law, absorption coefficients can be directly used as emissivity of dust grains in 
thermal emission spectra. 

The same technique has been employed to study the variations of band positions and band 
profiles with temperature, especially for low temperatures down to a few Kelvin. Because of 
the anharmonicity of the vibrational potentials of the lattice modes causing the infrared bands, 
the bands tend to get narrower and to shift to smaller wavelengths as the vibrational levels 
are de-populated with lowering T. Such studies have been performed e.g., by the groups in 
London, Kyoto, St. Louis and Jena for olivine particles 15214541 and for phyllosilicates (see 
Fig. [3]) 15514571 . It has been found that the sharpness and wavelength position of certain silicate 
far-infrared bands might be actually good indicators of the dust temperature. 

The drawback of the powder transmission spectroscopy in pellets, generally, is the possible 
influence of the polarisation of the embedding medium on the spectra. It can be easily shown 
for instance by Mie calculations that such an influence exists (see Sec. 3). In the general case of 
nonspherical or even irregular particle shapes, however, it is difficult to correct for this influence, 
although it is possible if the material's optical constants are known (see also below) 05811591 . 

4.2 Optical constants 

Most astrophysicists modeUng thermal dust emission spectra rely on calculated rather than ex- 
perimentally measured spectra, because of the flexibility of such calculations if parameters such 
as grain sizes have to be varied. The material properties enter into these calculations usually in 
the form of the so-called "optical constants". This term denotes the response function which 
characterizes the reaction of a material to an external (time-variable) electromagnetic field: the 
complex refractive index m ~ n + ik. If magnetic material properties can be neglected, m is 
the square root of the (also complex) dielectric function e which alternatively can be used as 
optical constants. 

It has to be noted that these quantities are macroscopic in character, therefore they loose their 
meaning for small clusters and molecules. In the transition region from solids to molecules the 



14 

12 

10 

M 8 
c 

6 
4 
2 

'io 30 50 100 200 

Wavelength (|a,m) 

Fig. 4. Far-infrared optical constants of calcite at temperatures of 10 K (solid lines) and 300 K 
(dashed lines) — after Ref . 6 1 . 

introduction of size-dependent optical material properties may be a reasonable way to include 
e.g., the confinement of charge carriers to the limited particle volume 1601 . 

The term optical constants is a bit misleading since they are strongly frequency-dependent. 
The dispersion with frequency is determined by resonances of the electronic system, of the 
ionic lattice and, at very low frequencies, by relaxation of permanent dipoles. In these frequency 
regions absorption becomes strong (high k) and n shows "anomalous dispersion", i.e. decreases 
with frequency. In many cases, this behaviour can be described by Lorentzian oscillators. 

Frequency-dependent optical "constants" of solid materials can be found in a number of 
databases which are available either in the form of books or electronic media. The most im- 
portant database in book-form is the "Handbook of Optical Constants of Solids" edited by E.D. 
Palik that currently consists of three volumes which appeared in 1985, 1991 and 1998 162^641 . 
These books are recommended since they comprise detailed discussions of the origin and errors 
of each dataset which may prevent the user from uncritical application of these data. 

The Jena-St.Petersburg-Heidelberg database of optical constants 1651 which is especially 
dedicated to cosmic dust applications can be found via the web pages ,http://www.astro.uni-jena.de| 
and |http : // www.mpia-hd.mpg . de/H JPDOC| 

As mentioned already, nearly all of these optical constants have been measured either on 
bulk samples or on thin films. The reason for that is that the planar geometry of films or 
coatings is much easier to describe in deriving optical constants from a measured spectrum than 
are irregularly shaped particles. However, even in the case of ideal bulk measurements, the 
determination of optical constants over a wide frequency or wavelength range is not a simple 
task. Since the material absorption in different spectral regions usually differs by many orders 
of magnitude, for the determination of k either transmission measurements on samples of very 
different thicknesses (from centimeter down to submicron scales) or transmission and reflection 
measurements have to be combined l66l . Many crystals show an anisotropy in their optical 
constants. In these cases, measurements with polarized light along the different axes of the 
crystal have to be carried out which require careful orientation of the crystal and alignment of 
the polarizers. 

In regions of strong absorption, as mainly considered here, usually reflection measurements 
have to be used and either the Kramers-Kronig relations (see Ref. 24) or a parametrization of 
the optical constants as a sum of Lorentzian oscillators have to be applied in order to obtain both 
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Fig. 5. Schematic structure of the aerosol spectroscopy setup. 



the real and imaginary parts of the optical constants. As already mentioned above, the optical 
constants of cosmic dust analogs may strongly depend on the temperature of the material, so 
that temperature-dependent data are desired. 

Despite of the difficulties inherent to the preparation of samples and the measurement and 
data evaluation procedures in obtaining such data, some progress has also been made in the 
measurements of the optical constants for dust materials within the last years. For instance, op- 
tical constants of amorphous magnesium silicates with a variable magnesium content have been 
measured in a wide wavelength range [67 1. For crystalline silicates, the optical constants of 
minerals of the olivine series (composition (Mg,Fe)2Si04) have been re-measured for different 
magnesium-to-iron ratios (e.g., Refs. 68 and 69) and even at temperatures down to 50 K (for 
Mg2Si04, ifTOl ). More temperature-dependent optical constants have been derived for carbon- 
ates ||6T1 , including the far-infrared wavelength range (see Fig.|4|i. 



4.3 Aerosol measurements and grain morphology 

The influence of the grain morphology and of the embedding into matrices has always been a 
concern in the use of experimental spectra for comparison with observed ones. As demonstrated 
in Sec. 3, the absorption band profiles of solid particles depend on the grain shape and on the 
electromagnetic polarisability (the refractive index) of the ambient medium. 

When comparing spectra of crystalline silicates measured in different laboratories, it has 
been realized that band profiles and peak positions showed systematic differences, which were 
supposedly related to the morphological properties of the particulates such as the grain shape 
(influenced by the powder grinding technology) and perhaps the (unknown) state of agglomera- 
tion inside the pellets 15411711 . On the other hand, spectra calculated from the material's optical 
constants assuming a spherical grain shape failed completely to reproduce both observed and 
measured spectra. Simple statistical approaches to more realistic grain shape models such as 
the Continuous Distribution of Ellipsoids (see Sec. 3) gave much better, but still not sufficient, 
results 1681 ■ 




Fig. 6. (a) Absorption spectra of irregular-shaped spinel particles measured in different media 
(nitrogen gas, KBr and CSi pellets - solid lines) compared to simulated spectra using the CDE 
(dotted line) and DFF (dashed line) models, (b) The shape factor distribution used for the DFF 
calculations (solid line) compared to that of the CDE (P(L) = 2(1-L), dashed line), (c) Electron 
micrograph of the particles. The black filter holes have a diameter of 0.5 yum. 



From the experimental side, an important step forward was made a few years ago, when 
the aerosol spectroscopy technique was introduced into the measurement of experimental com- 
parison spectra for astronomical puposes 117211 . In contrast to the pellet method, this technique 
proved to be suitable not only for avoiding any influence of an embedding medium onto the 
spectra, but also for a detailed study of morphological effects on infrared spectra of dust sam- 
ples 173)1 . The reason for this is that a direct electron microscopic investigation of the aerosol 
particles is routinely possible after the spectroscopic measurement. 

The experimental setup for aerosol infrared spectroscopy as used at the astrophysical labo- 
ratory of the AIU Jena consists of a White-type multi-path gas cell connected to a Bruker 11 3v 
FTIR spectrometer and equipped with an "aerosol generator" for dispersing fine powders in ni- 
trogen gas and filling the gas cell with this aerosol (see Fig.|5]i. The aerosol generator makes use 
of a Palas RBG 1000 rotating -brush particle disperser and is also equipped with an impactor- 
type grain-size discriminator in order to allow only small particles to reach the gas cell. The 
aerosol formed of these submicron-sized particles in a gas at about normal pressure is stable 
for several minutes, sufficient for FTIR measurements. After the measurement, a part of the 
particles are filtered from the aerosol for scanning electron microscopic (SEM) inspection. The 
drawbacks of this technique are that the column density of the particles is not precisely known, 
making a quantitative measurement of the absorption coefficient difficult, the large sample con- 
sumption, and the limitation to a certain particle size determined by the sedimentation speed. 

By comparing the aerosol-measured infrared spectra with pellet spectra, the influence of the 
embedding medium on the band profiles has been made clear for a large number of oxide fMl 
and silicate iTSl mineral "dust samples". The internet database making such aerosol spectra 
publically available (http://elbe.astro.uni-jena.de) contains also spectra of particles embedded 
in cesium iodide (Csl) pellets for comparison. The latter are usually characterized by infrared 
bands peaking at longer wavelengths. This does not mean that the bands are just shifted, but 
that surface modes closer to the transverse optical lattice resonance get a higher weight in the 
band profiles (see below). 

In the course of these aerosol measurements, great attention has been given to morphologi- 
cal effects in the spectra and it has been confirmed that grain shapes play a big role for the band 
profiles. For several compounds it has been possible to measure powders with different grain 
shapes, for instance more roundish ones vs. more irregular, sharp-edged shapes. The respective 
infrared bands were found to be characteristically different, for the roundish grains often reveal- 
ing a short-wavelength peak and a long-wavelength shoulder, whereas irregular grains produced 
long- or intermediate wavelength single peaks 05II74I . 

As mentioned in Sec. 3, the DFF model can successfully be used to reproduce these band 
profiles. Figure|6]demonstrates this for irregularly shaped spinel (MgAl204) particles measured 
in various environments (aerosol, pellets). The form factor distribution used in these computa- 
tions is also shown. It has to be noted again that this model, as many others, is only exact in the 
case of particles composed of an isotropic material. 

5 CONCLUSIONS 

The measurements and the theoretical calculations of optical properties of cosmic grains play 
an important role in a wide range of astrophysical applications. The collection of relevant 
data for well-characterized materials over a broad wavelength range has made great progress 
over the last two decades. Here we want to stress that the characterization of the physical and 
chemical structure of the materials is an important part of the study of cosmic dust analogs. An 
interesting direction of future research will be the comparison of the properties of cosmic dust 
grains collected by interplanetary space missions such as the Stardust mission (e.g., Ref. 76) 
with synthetic materials produced in the laboratory. 

The combination of condensation experiments carried out under conditions similar to cos- 
mic environments and advanced spectroscopic methods is certainly the future of solid-state 



astrophysics. Mass spectroscopy has to be used to characterize the transition from atoms and 
molecules to small clusters and solids. In addition, the change of material properties under the 
bombardment of cosmic rays and irradiation with intense X-ray and U V photons has to be better 
understood. Another important factor can be thermal annealing cycles of the grains. 

The application of laboratory data for the interpretation of astronomically measured fluxes 
rehes on the knowledge of the density and temperature structure of the considered dusty re- 
gions. Radiative transfer calculations have to be carried out for deeply embedded objects such 
as protostars and active galactic nuclei. In the case of very small particles the effects of quantum 
heating have to be taken into account. Even for the interpretation of spectra of optically thin 
configurations, the dust temperature and density distributions have to be known. 
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